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1 Introduction 



Motivated by the need to give answers to some fundamental questions in quantum field 
theory, during the last years there has been (and continues to be) a lot of interest in the 
problem of calculating the determinant of a differential operator, L (see for example Jl], 0). 
Often one has to deal in these situations with positive elliptic differential operators acting 
on sections of a vector bundle over a compact manifold. In such cases L has a discrete 
spectrum Ai < A2 < ... — 00. The determinant, detL = Hi Aj, is generally divergent and 
one needs to make sense out of it by means of some kind of analytic continuation. A most 
appropriate way of doing that is by using the zeta function regularization prescription 
introduced by Ray and Singer (see also 0, 0). In this procedure In detL is defined 
by analytically continuing the function A"'' In Aj in the exponent s, from the domain 
of the complex plane where the real part of s is large to the point s = 0. Introducing the 
zeta function associated with the spectrum Aj of L, 

i 

this is equivalent to defining 

IndetL = -C'(0). 

Only a few general methods for the exact evaluation of In det L are available. Thus, for 
example, given that the manifold has a boundary, in (see also [0, ||, |l^) the determi- 
nants of differential and difference operators have been related to the boundary values of 
solutions of the operators. When L is a conformally covariant differential operator, exact 
results may sometimes be obtained by transforming to a "more simple" operator L for 
which In det L is known. Then, the knowledge of the associated heat-kernel coefficients 
— nowadays available Hill , [l^ — gives sometimes the exact value of In det L [|l^, |14|, |15 



This approach has been used by Dowker to find the functional determinants for a variety 



of sectors of Euclidean space, spheres and flat balls for dimensions D < A |T6|, |17[ 
Similar techniques have proven to be very powerful in order to obtain estimates of different 
types 0, 11,11. 

As a rule, however, explicit knowledge of the eigenvalues Aj is necessary in order to 
obtain exact results for IndetL. This explicit knowledge of the eigenvalues is in general 
only guaranteed for highly symmetric regions of space, such as the torus, sphere or re- 
gions bounded by parallel planes. For these manifolds, detailed calculations have been 
performed in the context of Casimir energies and effective potential considerations (for a 



summary of results see [22 



In this paper we want to focus on a class of situations for which the eigenvalues of 
the operator are not known explicitly, but nevertheless the exact calculation of In det L 
is possible. The method developed is applicable whenever an implicit equation satisfied 
by the eigenvalues is known and some properties (later specified) of this equation are 
known. We exemplify our approach by taking L = —A on the /^-dimensional ball = 
{x G R'^; I a; I < R}, together with Dirichlet — or general Robin — boundary conditions. 
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In dimensions D < 5 a part of these results may also be obtained using the conformal 
transformation techniques mentioned above. (Here the restricition D < 5 results from the 
need of the knowledge of the heat-kernel coefficients a_D/2, which has been determined only 
recently for D = 5 [|12|.) This particular method was used by Dowker and Apps for D < A 



(this was before 05/2 were known) [^. Here we will apply a direct approach which has 
already been shown to be quite powerful for the calculation of the heat-kernel coefficients 
in the situations described above [|^. Interesting quantum field theoretical applications 
of the results obtained can be found in quantum cosmology p5[ [EB, HTl, E^, p9i . For these 



applications, the consideration of dimensions D higher than four are of interest, because 
the technicalities involved in dealing with higher-spin fields reduce essentially to the ones 
for scalar fields in those dimensions. Further applications are in statistical mechanics, in 



connection with finite size effects and in conformal field theory ||31|| . 

The paper is organized as follows. In section 2 we describe in detail our approach 
for the calculation of functional determinants using as an example the Laplace operator 
of the three-dimensional ball with Dirichlet boundary conditions. For generality and 
because the analytic continuation procedure employed is slightly easier, we start with the 
massive Laplacian, performing the limit m at a suitable point of our calculation. The 
result we find here agrees completely with the one recently given by Dowker and Apps . 
After having explained in detail the main ideas, we apply our approach to an (in principle) 
arbitrary dimension D and also to general Robin boundary conditions. We start in Sect. 3 
with Robin boundary conditions for the three-dimensional ball. The Neumann boundary 
conditions — a special case of the Robin boundary conditions — cannot be obtained as a 
limit of the parameter involved in Robin boundary conditions. An extra consideration, 
necessary in order to deal with this situation, is given at the end of section 3. In Sect. 4 
we describe in detail how our scheme can be applied in dimensions D > 3. Explicit results 
are given for dimensions D = 4,5,6 and all the different boundary conditions. The case 
D = 2 is briefly considered in Sect. 5. This case is slightly different because the lowest 
angular momentum / = needs to be treated in a speciflc way. In the conclusions (Sect. 
6) we summarize the main results of our investigation. As we will see in the course of 
our procedure, the functional determinant is naturally splitted up into two pieces. The 
contributions of each of these pieces and a detail of how they are obtained are given in 
three appendixes. 



2 Zeta function determinant on the 3-dimensional 
ball with Dirichlet boundary conditions 

In this chapter we want to concentrate on the 3-dimensional ball with Dirichlet boundary 
conditions in order to exemplify our direct approach, which is actually applicable in any 
dimension and for completely general Robin boundary condition. We are thus interested in 
obtaining the zeta function of the operator (— A + m^) on the ball = {x E R^ < R} 
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endowed with Dirichlet boundary conditions. The zeta function is formally defined as 



c(^) = EV, (2.1) 



k 



with the eigenvalues being determined through 

(-A + m2)0fc(x) = Afc0fc(a;) (2.2) 

{k is in general a multi-index here), together with the boundary condition. It is convenient 
to introduce a spherical coordinate basis, with r = |x| and angles = {6,(p). In these 
coordinates, a complete set of solutions of Eq. (|2.2| ) can be given in the form 



0z,m,n(r, fi) = r ■Ui^i{wi^nr)Yi^s{Q), (2.3) 



the J/+1/2 being Bessel functions and the 1^+3/2 hyperspherical harmonics [^. The wi^n 
(> 0) are determined through the boundary condition by 

Ji^{wi,nR) = 0. (2.4) 

In this notations, using A/^„ = wf^ + m^, the zeta function can be given in the form 

00 00 

c(^) = EE(2/ + i)Kn + ^')-^- (2.5) 

n=0 1=0 



Here the sum over n is extended over all possible solutions of Eq. ( |2.4| ) on the positive real 
axis, and (2/ + 1) is the number of independent harmonic polynomials in 3 dimensions. 

As it stands, the zeta function in Eq. is defined for 3?s > 3/2. The way how 

to construct the analytical continuation of Eq. ( p.5|) to the left has been explained in 
detail in Ref. [^| for the calculation of the heat-kernel coefficients associated with the 
operator (pl^). For this reason our description here will be brief. Starting point of the 
consideration is the representation of ({s) by the contour integral 

C(.) = E(2/ + 1) J ^ {e + mY'^lnJ,^.{kR). (2.6) 

^ — 'y 

where the contour 7 runs counterclockwise and must enclose all solutions of Eq. (|2.4| ) 
on the positive real axis (for a similar treatment as a contour integral see j^S], |3^, p5|). 
Subtracting and adding the leading asymptotic terms of Iv{zi') ioi v 00, v = I + 1/2, 
the following representation of C,{s) valid in the strip — l/2<3fJs<lis found to be valid 



(for details see p4 



2 



as) = 2Y^uZl{s)+Y. K{sl (2.7) 

1=0 i=-l 
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with the definitions 



sm vrs 



TT 



dz 



mRju 



In 



^^n[J.(.z.)] 



V2Tfu (1 



Dn{t) 



n=l 



(2.8) 



•2s oo 



E 



s+j 



Ch(2j + 2s-2;1/2), 



(2.9) 



and, for i = 1, 2, 



i2s oo 



Y: ^—jL^mRf^T{s + j)U2j + 2S-1; 1/2), (2.10) 



-fK E S^(mi?)^^CH(-l + ^ + 2j + 2.; 1/2) 
1-5 J }=0 3- 



i {i + 2a)T (^s + a + j + 
^0 ' r(l + a + |) 



Here 



Diit) 
D2{t) 



a=0 ° 

2 1 8^ S 



a=0 



t 16 



(2.11) 



and t = 1/vTT^, = a/1 + z"^ + ln(2;/l + x/TT^). Eq. ( p.7|) is a very suitable starting 
point for the calculation of the zeta function determinant C'(0); here we consider, for 
definiteness, m = 0. In the limit m only the j = term survives in the terms A^{s), 
n = —1, 0, 1, 2, and one immediately finds 



d_ 

ds 



E^ 



D, 



i=-l 



s=0 



3 ln2 Ini? 3 1 

~32 + ^ ~ ^ + 2^^ ^"^^ + 2^^ 



(2.12) 



For the part Z'^{s) some additional calculation is needed. First of all using the analyticity 
of Z'^{s) around s = 0, the derivative Z^'ifi) is found to be 



^d'(O) = - 



ln/,^(mz) — VT] + 



V 



(2.13) 



z=(mR)/v 
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and, in the limit m — 0, this reduces to 

1 1 

Z^'(O) = lnr(z/ + 1) + z/ - z/lnz/ - - \n{2nu) - —. (2.14) 

To perform afterwards the sum over z/, it is very convenient to use the integral represen- 
tation of lnr(z/ + 1) to find 

^£'(0)^/A(-^4-i + ;r^)^f. (2.15) 



Performing the sum, this yields 

Zo(0) = 2/* (-- + - + -^l^-^j— . (2.16) 

As is seen by studying the asymptotics for t — 0, the integral is well defined. For 
the explicit calculation of Z'^^^O), Eq. ( ^.Ifi] ), it is suitable to introduce a regularisation 
parameter and to define 



Z'^iO, z) = 2 I dtt^{-^ + -^+ ^-ri, (2.17) 



1 2 d e-' \ e-2 
2t2 + t^ + ^t(l-e-*)] 1^ 

with 

Z'^{0,0) = Z'^{0). (2.18) 
The individual pieces of the integral Eq. ( |2.17D may then be calculated by means of |^ 

and differentiating with respect to /5, 







oo 



r(i/ + i)^ / fx\ r(z/ + i)/i^ / ^ /i\ 



As a result we have 



Z'^iO, z) =Ch{z- 2; ^) r(2; - 2)[6 + z - z^] + \Ch [z; ^) r(^) 



and thus 



^D(0) = -^a-2)-^ln2. (2.21) 
Adding up the contributions from the Eq. ( |2.12| ) and ( |2.21| ) we end up with 

C'(0) = -| - ^1^2 - \c'^{-2) + \C'^{-1) - llni? (2.22) 
in agreement with Dowker and Apps [^. 
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3 Zeta function determinant on the 3-dimensional 
ball with Robin boundary conditions 

In order to treat Robin boundary conditions, only very little changes are necessary. Writ- 
ing the boundary condition in the form 



— Jl+l/2iyJl,nR) + Wi^riJUl/2i^l,nr) \r=R = 0, 

the starting point of the calculation, analogous to Eqs. (|2.8|)-(PTTT|), is 



Sm TTS 



TT 



oo 
















R/v 







-In 



n=l 



V" 



(3.1) 



(3.2) 



D, 



and for i = 1,2, 



T32s oo 

■fTT E ^—^{mRf'CH{-l + 1 + 23 + 2s; 1/2) 
2^ (z + a)r(g+j + ^ 



a=0 



r(i + ^ 



Here we need the polynomials 



and 



M2{t,u) = ^2;2,at 



2+a 



a=0 



— + U+ U — i 

16 2 2 / V8 2/ 16 



The contribution from the asymptotic terms is 

d 2 



ds ^ * 



Ri 



i=-\ 



s=0 



1 ln2 Ini? 3 1 

32 ~ ^ + ^ + 2^^ ^"^^ " 2^^ ^"^^ 



+- +-iu^ + 2u^ ln2 + u^ Ini?. 
2 



(3.3) 
(3.4) 



(3.5) 



(3.6) 



(3.7) 



(3.8) 
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For the rest, things are very similar to the Dirichlet case. We find 

oo 

Zn'iO) = Idt [--+^- -- + -—]— - (3.9) 





1 


1 


1 




' 2 ~ 


1 




" ) 


u 








+ - 

V 




2z/2' 



+ ln 

and thus, using Eq. ( |2.15| ), we have 

Z'M = -^C«(-2) - ^ ln2 + 2f j-z/ln (l + ^) + n - . (3.10) 

The remaining sum may be done by differentiating with respect to u and integrating back. 
As a result 

(3.n, 

= -2"^(7 + 21n2) -MlnT (2 + + / (2 + ^) • 



Let us mention that the last equation can also be given in terms of the Barnes G-function 
Putting things together, we arrive at our final result, 

C(0,u) = —--\n2--Cf,(-2)--CR(-l) + —\nR (3.12) 

u 

+^ - 2Mlnr (2 + w^lni? + 2 J dx InT (2 + ^) ' 



which concludes the consideration of the three dimensional case with Robin boundary 
conditions. The dependence of C'(0, u) on the parameter u is shown in Fig. 1 for R = 1. 

Eq. ( |3.12| ) can be given in an alternative form using the following expression, first 
obtained in 13611 (see also 1221) 



j dx \tit(^- + x^ = Ch{-1,u + 1/2) + Ch'{-1,u + 1/2) 


+ |ln(2vr) + i^^ + iC/?'(-l)- (3.13) 

As is seen, the limit u — >■ —1/2, corresponding to Neumann boundary conditions, is 
not smooth, since a logarithmic divergence appears due to the lnr(l/2 + -u) term. This 
logarithmic divergence might be traced back to the v = 1/2 term in Eq. ( |3.9|) . Thus, in 
order to find the functional determinant corresponding to Neumann boundary condition. 
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zeta' (0, u) 




Figure 1: Plot of the dependence of C'(0, u) on the parameter u for R = 1 and 
for dimensions D = 3,4,5,6. Notice the divergence that apppears 
for u = 1 — D /2 in each dimension, corresponding to the case of 
Neumann boundary conditions. 



the term = 1/2 has to be treated separately. In fact, looking in detail at the derivation 
of Eq. ( p. 91) starting from Eq. (|3.2|) , it is seen, that the case u = —1/2 has to be treated 
specifically, because the behavior of {u/ R)Iy{zi') + {zv / R)Il{zv) for 2; — >• is different for 
u = -1/2. 

Probably the easiest way to find the results for the Neumann boundary conditions is 
to write 



C(s,-l/2) = C='{s, -1/2)+ hm (C{s,u)-C='{s,u 

M— >— 1/2 ^ 



(3.14) 



because then we can use all the results for the Robin boundary conditions that we have 
derived before. Here C'^°(s, u) is the contribution from the angular momentum component 
/ = to u), 



c 



1=0, 



s,u 



Sm TTS 



TT 



/ dk -m^l-^^ln 
J ok 



-h/2{kR) + kI'y,{kR) 



(3.15) 



Proceeding with the calculations as for the Robin boundary conditions, one easily finds 
d 



ds 



(C'=°(.,-l/2)-C'=°(5,n)) 



s=0 



-21ni2 + ln(3/2) + ln(u + 1/2). 



(3.16) 
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In the limit m — — 1/2, the logarithmic divergence in Eq. ( p.l6| ) cancels the divergence in 
({s, u), Eq. (|3.12| ). Thus the limit m — > —1/2 is well defined and the functional determinant 
for Neumann boundary conditions reads 

m-i/2) = -l-lM + i„3-«'"« 3««'<-2) 



32 6 24 4 2 

1/2 

-2 / dx lnT{x). (3.17) 







This concludes the consideration of the 3-dimensional case: all ordinary boundary condi- 
tions have been dealt with explicitly. As for the Robin boundary condition, an alternative 
form may be presented using 

1/2 

/ dx lnr(x) = - + — In2 + -In7r--Cij'(-1). 
J ^ ^ 8 24 4 ^ 



4 Zeta function determinants on D > 3 dimensional 
ball 

For higher dimensional balls, exactly the same procedure may be employed in order to 
calculate the zeta function determinant. The starting point is here (for details see ||24]| ) 

°° r dk f) 

Us) = T.MD) J — {e + mY'^lnJ,^n^{kR), (4.1) 

ry 

with the number di{D) of independent harmonic polynomials, explicitly 

di{D) = {2l + D-2) ^\^^^-^f . (4.2) 

It is suitable to introduce the coefficients ea{D) by 

D-2 . D-2\" 

di{D) = ea{D) (l + . (4.3) 

a=l ^ ^ / 

A representation which shows the analytic structure around s = may then be given in 
analogy with Eq. ( |2.7| ) (now one defines u = I + [D — 2)/2) 

oo D-1 

as)=Y.di{D)ZUs)+Y: A^is). (4.4) 

/=0 i=-l 
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Here we have used the following definitions. First 



Zl(s) 



Sm TTS 



TT 



dz 



mR/u 



ZV 

'r 



m 



d_ 

dz 



[Inl^^uz) — vrj 



D-l 



Dn{t) 



+ \n{V2^{l + z^fl^)-Y.^ 



n=l 



(4.5) 



The polynomials Dn{t) arise from the asymptotic expansion of Iv{vz) In detail one 
has 



iJvz) ~ 



1 



1 + E 



Uk{t) 



k=l 



(4.6) 



with the recursion relation 



Uk+i{t) = -t'{l - r)4(t) + - / rfr (1 - 5r^)M,,(r) 



(4.7) 



for the polynomials Uk(t). The coefficient functions Dn{t) are then defined through the 
cumulant expansion 



In 



k=l 



^ Dnjt) 
,,n 



(4i 



n=l 



and are easily found with the help of a simple computer program. 

The Af{s) have already been determined in and we give only their final form for 
completeness. Introducing 



they read 



a=0 



i+2a 



(4.9) 



d2s oo / ly- Y (j + s - ^ 



4v^r(s)^o j! 



s + j 



-D-2 

J2 eaiD)CHi2j + 2s-l-a;{D- 2)/2) 

.a=l 

j^2s oo 



4r(s)§ j! 



1)^ 



i?)^^r(s + j) 



D-2 



J2 ea{D)CH{2j + 2s-a;{D- 2)/2) 



.a=l 



(4.10) 



(4.11) 
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R ^ 



-iV 



2T(s 



niR) 



2j 



j=o J- 

E ea{D)CHi-a + t + 2j + 2s; {D - 2)/2) 



(4.12) 



.a=l 



X E/ "^iii" 
a=0 



+ 2a)r(s + a + j + 



r(i 



a + 



Their contributions to the functional determinant are easily determined and are listed in 
Appendix A, for the dimensions D = 4,5, 6. 

For the calculation of Z'^[s) we go on as for the three dimensional case. We have 
found 



1 



D-l 



Z^'(O) = lnr(z/ + l) + z/-z/lnz/--ln(27rz/)+ E ^ 



n=l 



1 1 1 





°° /D-l 



1/ t 



+ E 

n=l 



tAn-l)\ ^2 t^d-l) t 



tu 



(4.13) 



This is the suitable starting point for the summation of the angular momentum. Intro- 
ducing the regularized version of Z£){s) as in Eq. ( ^.17] ) one obtains 



D-2 



— d'^ n 1 



D-2 



e-* dt" V t e* - 1 



+ Ee.(^) E 



DM) 



a=l 



n>Q+l 



T(n — a) 



Ch [ z + n - a; 



D-2 



T{z + n — a) 



D-2 



-iria + iy.CHiz-a-l; 



r{z — a 
D-2 



(4.14) 



rf^; — a — V 



which we need at z = 0. The remaining task is the calculation of the integral term in 
Eq. (4.14). It may be given by repeated differentiation of Eq. (2.19). Defining 



1 g '^x 

1 - e-^^ 



dx 



and 



f{k,p,-f,u) 



X' 



u+k-2 -[^+{k-l)l3]x 



[1 - e-P^Y 



(4.15) 



(4.16) 
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we get the recurrence 



d 



k-ldp' 



(4.17) 



The first few cases, k = 2,3, 4, and a general formula for f{k, (3, 7, v) is given in Appendix 
B. All terms appearing in Eq. ( fl.l4| ) and resulting from the integration may be given in 
terms of this functions f{k,(3,'^, v). In detail this is seen as follows. First of all one may 
show, that 



1 1 



dt" V t e* - 1 
with the recursion relations 



k+l 



(k) 



ik—l—a 



k=0 



1=1 



ly 



(a) 



-ira\, 



B 



(a) 



{k - a)Bt'^ + Bf- 



?("-!) 



for /c = 1, a — 1, 



for B^f* , and similar ones for AY' H 



(k) 



A 



(k) 

k+l 



(fc) 



(fc-i) 



for 2 <i<k, 



coming from 



k+l 



A 



(k) 



\ (e* - 1)^ 



As a result, 



D—2 °° 

^ e^{D) J dt t^^ 

a=l n 



e-*rft" Vt e* - L 

D-2 a k+l , n _ 9 



a=l 



k=0 



i=l 



z + k — a — i] . 



Thus we have derived all necessary equations, showing that Z'j^{0, z) may be given solely 
in terms of F-functions and Riemann zeta functions and the limit 2; — may be taken. 
The results are listed in Appendix C for D = 4, 5, 6. 
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Adding up the contributions from the asymptotic terms and Z'^{0), we have found 
the following final results for the zeta function determinants of the Laplace operator with 
Dirichlet boundary conditions in D = 4, 5 and 6 dimensions, 



cm 



CM 



a(o) 



1 T'? 1 1 1 

30240 90 90 3^^^ ^ 

-^a-2) + ^C«(-l), 



47 17 , , 



1 7 ^ 
5760 64 ' 



9216 2880 

+— Cr(-3) - — Cr(-2) - — Cr(-1) 
48 ^ ; 32 ^ 48 ^ 

4027 1 . 1 . „ 1 



6486480 756 756 



In 2 + — In i? + - C;j(-5) - - C'r{-^) 



60 



24 



24 ' 60 

The result for D = A agrees with Dowker and Apps ||2^. As is clear from our presenta- 
tion, any higher dimension D can be treated in exactly the same way without additional 
problems (they just get a bit arithmetically cumbersome). 

Let us now describe briefly the calculation for Robin boundary condition. There the 
starting point is 



as,u) = Y.di{D) J — (k' + mT— In ^-J,(kR) + kJUkR) 



(4.18) 



The relevant asymptotic expansions are ( [4. 6] ) together with 



€^^^(1 + ^2)1/4 



1 + E 



Vkit) 



' ' ' \J2txv z 
The relevant polynomials this time are determined by 



fc=i 



In 



k=\ 



k=\ 



E 

n=l 



M„(t) 



where, in analogy to Eq. 



2i 



), we define 
M„(t,M) = 5]z,,,f+». 

a=0 

Continuing as for the 3-dimensional case, and in analogy to Eqs. 
sm vrs 



TT 



dz 



mR/ V 



\r) 



In 



m 



oz 



u 



(4.19) 
(4.20) 

(4.21) 
( [4.5| ), we define 
(4.22) 



zv 



RV2 



l + z')-^ 
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In the limit m — 0, we obtain 



r /II 1 \ e^''" 
Zh'W = /* + ^ — (4.23) 





and, realizing that 



M„(l, 0) = M„(l, m) - M^l, 0) = (-1)"+' — , (4.24) 

n 

we find, after performing the sum over the angular momentum, 

Z'^{Q) = Z^(0) (4.25) 

+ E E e.(^)^" |- In (i + ^) + E (-1)"^'- (-)l • 

The remaining sum may once more be done by differentiating and integrating back, and 
we end up with 

Z'^{Q) = Z^(0) (4.26) 

f /D-2\ ?y"+i 



2 / a + 1 
+ E {-ir^\H{n-a;{D-2)/2)- 

n=a+2 ^ 
1- /D - 2 

- (-l)"a y dx x"-^ In T f + x 



+(-ir«"inr(^ + «)}. 

The results for the special dimensions D = 4,5 and 6 are not listed, because they are 
essentially given in Eq. (|4.26|) together with Appendix C. Let us only note that 

ei(4) = 0, 62(4) = 1, 

ei(5) = 62(5) = 0, 63(5) = ^, 

ei(6) = 0, 62(6) = -^, 63(6) = 0, 64(6) = ^. 

The contributions of the asymptotic terms are calculated from 

A^,{s)=A^,{s), A^{s) = -A^{s), (4.27) 
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and 

D-2 

^ ea{D)CHi-l + 1 + 2J + 2s; (D - 2)/2) 

a=l 

2i (i + a)r(s + j + ^^) 

They are collected in Appendix A. 

Summing up, we have found the following final results for the functional determinant 
with Robin boundary conditions: 

11 u 5u^ ln(2) ;/;Mn(2) 

~ 4320 ^ 30 " T2" ~ y ^ "90" ^ 3 

ln(i?) InjR) Cr{-^) Cl^(-2) (^^(-1) 

90 3 3 2 6 

u 

Wlnr(l + 'u) -2 J dx xlnr(l + x), 



(5(0, 





61 11m 11 



46080 576 16 72 24 

7 ln(2) 17 In(i^) ln(i?) ln{R) 

^ 720 5760 24 ^ 12 

5C^.(-4) 7C;,(-3) Cr{-2) ^C'ni-^) 



64 48 32 48 

X 



J dx In r f - + a; j + J dx x? In F 



12 



+— It In r f- + - -M^ In r f - + M 

12 V2 ; 3 V2 
, _ _ 9479 u hVIu^ 83 19 



32432400 315 15120 1512 480 

log(2) ln(2) ^5 ln(2) ln(i?) In(i^) 

'45 756 36 ^ 60 ^ 756 ^ 36 

\^{R) , C;,(-5) , Ck(-4) C;,(-2) Clil-l) 



+ 



+ 



60 60 24 24 60 

1 " 1 " 

-^^ Ax xlnr(2 + x) -:^Jdx x^lnr(2 + x) 



-^u'^\nr{2 + u) + —u^\nr{2 + u). 

The detailed dependence of ('{0, u) on the parameter u for dimensions L> = 4, 5, 6 is given 
in Fig. 1, for R^l. 



16 



Finally, we are left with the task of the calculation of the zeta function determinant for 
Neumann boundary conditions. As we have already seen, the / = term corresponding 
io u = [D — 2)/2 has to be treated separately. Employing the same procedure as for 3 
dimensions, one finds 



ls=0 



Y.^o.{D) (^^y [-2\nR + \n{D/2) + ln((D - 2)/2 + u)] 

which results in the following final results for Neumann boundary conditions, 

493 61 1n2 151 In i? 

(^'(0, -1) = + 

^ 4320 90 90 

I a(-3) ^ C'Ri-2) , C'Bi-^) 



3 2 6 

1 

+2 J dx {x- l)lnr(a;). 



w,^ o/oN 19261 713 In 2 , ^ 9647 In i? 

(5 (0,-3/2 = h n5 

^ ' ' ^ 46080 720 5760 

7Ck'(-3) C/?'(-2) ^ Ci?'(-1) 



64 48 32 48 

3/2 3/2 

X — - j Inr(x) 



+— j dx Inr(x) — j dx 



7087979 1181 ln2 , 6379 In i? 

Ce 0, -2 = + ln3 

' ^ 32432400 3780 3780 

Cr'{-5) Cn'{-4) _ a{-2) _ 

60 24 24 60 

2 2 



-- / dx (x-2)lnr(x) + - / dx (x - 2)^ In r(x). 
6 J 3 J 



This terminates the explicit calculation for all kind of boundary conditions and for dimen- 
sions D = 3,4,5,6 (that can be extended immediately to any dimension D > 6). Once 
more the results may be given in an alternative form by using Eq. (|3.13|) , but this makes 
little difference for numerical calculations. 
Finally, we are left with the case D = 2. 
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5 Zeta function determinants on the 2-dimensional 
ball 



For the 2-dimensional ball the procedure has to be changed slightly. Here the degeneracy 
of every Z > 1 is 2, / = has to be counted only once. Due to the presence of this term 
I = 0, the starting point Eq. (|2.7| )- (|2.11 ) is not valid any more and may be applied only 
to / > 1. The / = term may be treated as before for the Neumann boundary conditions. 
We shall not give any further details for this case but only write down the final results 



— which on the other hand have appeared in part in the literature |38|. We quote them 
for completeness. 

For Dirichlet boundary conditions we have 

CD(0) = ^ + 2a-l) + ^ln7r + ^ln2 + ^lni?. 

The zeta function determinant for general Robin boundary conditions reads 

Cfl(O) = --^ + ^lni? + 2CRi-l) - ^ In 2 - ^ In TT + 2m ln(2/i?) - In m + 2 In 1(1 + m). 

Finally, the results for Neumann boundary condition is 

G(0) = - ^ Ini? + 2C;,(-1) + ^ ln2 - 1 Invr. 

And this concludes the list of examples of zeta function determinants on the ball that we 
had promised to consider. 



6 Conclusions 

In this paper we have developed a systematic approach for the calculation of functional 
determinants of (elliptic differential) operators, which is very useful in all cases when 
the basis of functions — constrained by the equations cooresponding to the boundary 
conditions — is known. Using our approach we have calculated the zeta function determi- 
nant of the Laplacian on the ball for Dirichlet, Neumann and the general Robin boundary 
conditions. Explicit results for dimensions D < 6 have been given. All necessary formulas 
to iterate the calculation of C'(0) ^"^^ obtain it in any dimension, by means of a simple 
computer program, have been given explicitly. 

An extension of the present work to higher- dimensional bundles, such as those for 
spinors and vectors, is envisaged. This should allow to establish connections with recent 
work dealing with mixed boundary conditions B^, HO, Hi . 



Comment: Using a different approach, based on work by Moss |^ and Voros 
Dowker has also considered the calculation of heat-kernel coefficients and functional de- 
terminants for the Laplace operator on the ball with Robin boundary conditions |44]. We 
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are indebted with him for interesting and fruitful correspondence during the course of the 
present work. 
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A Appendix: Contributions to the zeta function de- 
terminant from the asymptotic terms 

In this appendix we list the contributions of the asymptotic terms to the zeta function 
determinant for Dirichlet and Robin boundary conditions. 



A.l Dirichlet boundary condition 

For Dirichlet boundary conditions we find 



83 7 , ln2 InR ^, , Cr( 



{Asym)D[D = ^) = Cr(-3) - 

^ ^ > 6048 360 90 90 ' 

47 TT^ 11 17 

(Asym)D(-D = 5) = + In 2 + Ini? 

^ ^ > 9216 8640 5760 5760 



+^c;^(-4) + ^c^(-3) - lc;,(-2) - ^c;,(-i), 

40091 7 In 2 Ini? Ck(3) Cr(-5) 



25945920 3360 756 756 15120 12 

c;,(-4) ^ c;,(-3) ^ c;^(-2) 



24 12 24 

A. 2 Robin boundary condition 

For Robin boundary conditions the final results read 

73 1 1 1 1 

{Asym)R{.D = ^) = 7+— In 2 Ini? - CU-3) + -Cr(-2) 

\ y )R\ ) ^32Q 360' 90 90 ) ^ ^^R\ ) 
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(Asym)„(D = 6) 



111 1 1 1 1 

u H u 7M H u H vr M'^ H u 

576 48 24 ^ 24 72 24 



12 ' 12 6 



243079 1 



— In 2 ^Ci?(3) 

756 15120^^ ^ 



129729600 3360 



12 
1 



24 



12 



24 



315 



u + 



517 
15120 



^24 ^5 
TT M U 

288 45 



41 

1512 



1 



36 60 60 ^ 



B Appendix: Detail of the calculation of Z'jj{0) 



Using the recurrence relation ( ^4.17| ), the function f{k, j3, 7, v) for = 2, 3 and 4 read, 

roo 3;i'g-(7+/3)a: 



/(2,/5,7,i^) 



(1 - e-/5^)= 



r(^ + i) 



/(3,/3,7,i^) = 

_ r(z/ + 2) 
2/3^^+2 

/(4,/5,7,i^) = 

_ r(z/ + 3) 
2(3^+^ 



^i/+lg-(7+2/3)x 



(1 - e-^^)3 



0/ ^ 'y 'y 'y 'y 

C(..^)-(i + 22)c(. + i.I) + 2(i + 2k(. + 2.1) 



00 2;'^+2g-(7+3/3)x 



(1-e- 



■ da; 



C(z/, I) - 3(1 + 2) C(z/ + 1, ^) + [2 + 6^ + 3(^)2] + 2, ^ 



/5' 



For the general formula one has 

T{v + k-l) 



where 



7 

Cfc-i,o(-^) 



7^ 



1, c._.,(^) = ^ + (^ + l) + -- + (^ + fc-2), 



(2.1) 



(2.2) 



(2.3) 



(2.4) 
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,7, _ + 

C Appendix: Contribution of Zd{s) to the zeta func- 
tion determinant 

In this appendix we list the contributions of Zd{s) to the zeta function determinant for 
dimensions D = 4,5 and 6. They are 
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Figure caption: 

Figure 1: 

Plot of the dependence of ^'(0, u) on the parameter u ior R — 1 and for dimensions 
D — 3,4:, 5, 6. Notice the divergence that apppears ior u — 1 — D/2 in each dimension, 
corresponding to the case of Neumann boundary conditions. 
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